Existence of Chaos for a Singularly Perturbed NLS Equation 

Yanguang (Charles) Li 

Abstract. The work ]]J is generalized to the singularly perturbed nonlinear 
Schrodinger (NLS) equation of which the regularly perturbed NLS studied 
in [Q is a mollification. Specifically, the existence of Smale horseshoes and 
Bernoulli shift dynamics is established in a neighborhood of a symmetric pair of 
Silnikov homoclinic orbits under certain generic conditions, and the existence 
of the symmetric pair of Silnikov homoclinic orbits has been proved in 
The main difficulty in the current horseshoe construction is introduced by the 
singular perturbation ed% which turns the unperturbed reversible system into 
an irreversible system. It turns out that the equivariant smooth linearization 
can still be achieved, and the Conley-Moser conditions can still be realized. 



1. Introduction 

Consider the singularly perturbed nonlinear Schrodinger (NLS) equation, 

(1.1) iq t = q C( + 2[M 2 - u 2 }q + ie[q (( -aq + fi], 

where q = q(t, £) is a complex-valued function of the two real variables t and £, 
t represents time, and £ represents space. q{t, £) is subject to periodic boundary 
condition of period 2-7T, and even constraint, i.e. 

q(t, C + 2tt) = q{t, Q , <?(*, ~C) =<?(*, • 
a > and (3 > are constants, and e > is the perturbation parameter. For 
simplicity of presentation, we restrict w by u> G (1/2,1). In Q, the following 
theorem on the existence of Silnikov homoclinic orbits was proved. 

Theorem 1.1. There exists a cq > 0, such that for any e € (0,eo), there exists 
a codimension 1 surface in the external parameter space (a, /?, u) € M + x M + x 
R + where ui G (i,l)/S, S is a finite subset, and aoj < (3. For any (a,f3,uj) on 
the codimension 1 surface, the singularly perturbed nonlinear Schrodinger equation 
( |l.l| ) possesses a symmetric pair of Silnikov homoclinic orbits asymptotic to a saddle 
Q e . The codimension 1 surface has the approximate representation given by a = 
1/k(u>), where k(u>) is plotted in Figure^. 

Notice that if q(t, £) is a homoclinic orbit, then q(t, ( + ir) is another homoclinic 
orbit. Thus q(t, £) and q(t, C + tt) form a symmetric pair of homoclinic orbits. Based 
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Figure 1. The graph of k(u>). 



upon the above theorem, we will construct Smale horseshoes in a neighborhood of 
the symmetric pair of homoclinic orbits. The construction is a generalization of 
that in JIJ where the singular perturbation ed%. is mollified into a bounded Fourier 
multiplier. The main difficulty in the current horseshoe construction is introduced 
by the singular perturbation edfj, which turns the unperturbed reversible system 
into an irreversible system. Specifically, denote by F* the evolution operator of 
the singularly perturbed nonlinear Schrodinger equation ( |l.l| ). When e = 0, Fq is 
a group. When e > 0, F* is only a semigroup. It turns out that the equi variant 
smooth linearization can still be achieved, and the Conley-Moser conditions can 
still be realized. Of course, one has to replace the inverse of the evolution operator 
.F* by preimage. The article is organized as follows: In section 2, we present 
equivariant smooth linearization. In section 3, we present the Poincare map and 
its representation. In section 4, the fixed points of the Poincare map is studied. 
In section 5, we present the existence of chaos. Finally, in section 6, numerical 
evidence for the generic conditions is presented. 

2. Equivariant Smooth Linearization 

The symmetric pair of Silnikov homoclinic orbits is asymptotic to the saddle 
Q € — ^/le %e , where 

(2.1) / = W 2 -ei#^V + -, CO s0 = ^, 0€(O,~). 
Its eigenvalues are 



(2.2) xt = -e[a + n 2 ]±2\l{-+^-I){ZI-^--) , 



n 2 



where n = 0, 1, 2, • • • , co 6 (i, 1), and I is given in (2.1). The crucial points to 
notice are: (1). only Aq~ and Xf have positive real parts, RejAp } < Re{A^}; (2). 
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all the other eigenvalues have negative real parts among which the absolute value 
of Re{A^} = Re{A^} is the smallest; (3). |Re{A^}| < Re{A^}. Actually, items (2) 
and (3) are the main characteristics of Silnikov homoclinic orbits. 

Lemma 2.1. For any fixed e € (0, eo), let E e be the codimension 1 surface in the 
external parameter space, on which the symmetric pair of Silnikov homoclinic orbits 
are supported (cf: Theorem 1.1). For almost every (a,/3,w) £ E e , the eigenvalues 
A^ satisfy the nonresonance condition of Siegel type: There exists a natural 

number s such that for any integer n > 2, 



> 1/r 



3 = 1 

for all r = 2, 3, • • • , n and all l±, 1%, ■ ■ ■ ,l r G Z, where A„ = A^ for n > 0, and 
A„ = \Z n -i for n < 0. 

Proof. The same proof as in pi can be carried through here. QED 



Thus, in a neighborhood of Q ei the singularly perturbed NLS (1.1) is analyt- 
ically equivalent to its linearization at Q e ||. In terms of eigenvector basis, ( |l . l[ ) 
can be rewritten as 

x = -ax - by + j\f x (Q), 
y = bx - ay + j\f y (Q) , 
(2.3) h = 71-21 +A4 1 (<9), 

Q - LQ+jV q (Q); 

where a = — Rc{A^}, b = Im{A^}, 71 = A^, 72 = A^ ; Af's vanish identically in 
a neighborhood Q of Q — 0, Q — (x,y,Zi,Z2,Q), Q is associated with the rest of 
eigenvalues, L is given as 

LQ = -tQ Ci: - 2i[{2\Q e \ 2 - uj 2 )Q + Q 2 Q] + e[-aQ + Q cc ] , 

and Q e is given in ( |2.l| ) . The following theorem on well-posedness is standard Q . 
Let F* (0 < t < 00) be the evolution operator of the singularly perturbed NLS 
(P^), and H s be the Sobolev space. 

Theorem 2.2. For any s > 1, and any Q E H s+2 , F*(g ) € C°([0, 00); H s+2 )n 
C 1 ([0, 00); H s ). For any fixed t € [0, 00), F t is a C 2 map in H s . 

3. The Poincare Map and Its Representation 

Denote by hu {k = 1, 2) the symmetric pair of Silnikov homoclinic orbits. The 
symmetry a of half spatial period shifting has the new representation in terms of 
the new coordinates 

(3.1) a o (x, y, zi,z 2 , Q) = (x, y, Zi,-z 2 , a o Q). 

We have the following facts about the homoclinic orbits: 

1. The homoclinic orbits are classical solutions, 

2. As t — * —00, the homoclinic orbits are tangent to the positive zi-axis at 
Q = 0. 

The same proof as in Q works here for item (2). Since a is the smallest attracting 
rate, we assume that 
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• (Al). As t — > +00, the homoclinic orbits are tangent to the (a;, y)-plane at 
Q = 0. 

The Poincare section is defined as in Q . 

Definition 3.1. The Poincare section So is defined by the constraints: 

y = 0, ryexp{— 2na/b} < x < 77; 

< z\ < 77, -rj < z 2 <r], \\Q\\ < 77; 

where 77 is a small parameter. 

The horseshoes are going to be constructed on this Poincare section. The 
auxiliary Poincare section is defined differently from that in (5J. 

Definition 3.2. The Poincare section Si is defined by the constraints: 

Z\ =77, -77 < z 2 < ?/, 



^x 2 + V 2 < V, ||Q||<?7- 
The Poincare map is defined as follows. 
Definition 3.3. The Poincare map P is defined as: 

P : [7cS ^S , P = P 1 °°P 1 , 

where 

Po 1 : C/oCEowEi, VQeC/o, P X (Q) = e Si, 

and to = to(Q) > is the smallest time t such that F t (Q) £ Si, and 

P? : C/ 1 cSi^^(=SoUaS ), VQGC/i, e^>, 
and ti — ti(Q) > is the smallest time t such that F*(Q) £ So. 



The map Pq 1 has the explicit representation: Let Q° and Q 1 be the coordinates 
on S and Si respectively, y° = 0, and z\ = 77, then to = — In and 



a; cos 



2; sin 



b 77 

— In — 

& 77" 

— In — 

7i *i 



Q 1 



= e toL Q° 



Let and be the intersection points of the homoclinic orbit hi with So and 
Si respectively. The discussion with respect to the other homoclinic orbit h 2 is the 
same. In a small neighborhood of Q\, the map P® has an approximate representa- 
tion. By virtue of the fact that the homoclinic orbit hi is a classical solution, the 



Well-Posedness Theorem 2.2 implies that F t (Q\) is C 1 in t. Thus for Q 1 in a small 
neighborhood of Ql, 

(3.2) P^Q 1 ) = P°(Ql) + CiQ 1 - Ql) + OiWQ 1 - Qlf), 
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where 
CiQ 1 - 

and t\ = 
vanishes, 

Thus 



(3.3) 
Let Q 

(3.4) 

where 



Q\) = d^{Q\) o (Q 1 - Q\) + d t F^(Ql) o -±-(Qt) o (Q 1 - Q\) , 
tiiQ 1 ) is defined by the constraint that the y-coordinate of F tl {Q 1 ) 



dQF^{Q\) + d t F^{Q\)-^{Ql) = Q 
dQ 1 



Q° — , and Q = Q 1 — Q\, then Pf has the approximate representation 



/ 5° 

?° 


\ 




( x 1 

v 1 


\ 




= c 




~A 






~4 






J 




K Q 1 


) 



Z-O (xr + (yr + &y + \\Q 



( Cn 


Cl2 


Cl3 


C14 \ 


C21 


C22 


C23 


C24 


C31 


C32 


C33 


C34 


V c 41 


C42 


C43 


C44 y 



in which Cj; (j, 1 = 1,2, 3) are real constants, Cj4 (j — 1, 2, 3, 4) and C41 (I = 1, 2, 3) 
are linear operators. 



(4.1) 



4. The Fixed Points of the Poincare Map P 

As to — > +00, to the leading order, the fixed points of P satisfy 

/ x° \ ( xlcosbto \ 





c 



x® sin bto 
V J 



where 



*1 — p atojl a0 _ ai o ~0 /)0 _ at /SO 
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Explicitly, the second and the third equations in ( |4. l[ ) are: 
%1 c 2 i cos bt Q + c 2 2 sin bt Q 



(4.2) 



c 3 i cos bt + c 32 sin &£ c 



C23Z\ = 0, 
- c 33 zl = 0. 



Lemma 4.1. C23 and C33 do noi vanish simultaneously. 

Proof. Notice that W u (Q e ) is two-dimensional, and intersects So (or its exten- 
sion to — 77 < z\ < rj) into a one-dimensional curve with tangent vector 

( °\ 


1 

V / 



v =C 



Notice also that for any Q G hi, 

(4.3) dim{T W u (Q e ) n 7^W s (Q e )} = 1, 

where 7^ denotes the tangent space at Q. If C23 and C33 vanish simultaneously, 
then v £ T^ W s (Q e ) which implies that 

dim{T go W u (Q e ) r\T d „W s (Q 6 )} = 2 



which contradicts (4.3). The lemma is proved. QED 
Let 

Al = C21C33 - C31C23, A 2 = C22C33 - C 32 C23- 

We assume that 

• (A2). Ai and A2 do not vanish simultaneously. 



Then ( |4.2| ) has infinitely many solutions: 



(4.4) 
where 



4' 3 = ^[lv-<pi], 



1 



ipi = arctan{Ai/A2}. 



Without loss of generality, we assume C23 7^ 0. Then, solving Eqs. (|4.2| ), we have 
(4.5) zf' l) = -^[C2 3 ]" 1 {c2icos6i ° +c 2 2sin6i °}. 



Solving (4.1), we have 



(4.6) 
(4.7) 



Q 



(0,1) 



en cos Wq + C12 sin bt\p 



C41 cos frig + C42 sin bt^ 



Cl3^ 2 1 



+ ^43 ^2 



Finally, by the implicit function theorem, there exist infinitely many fixed points of 
P, which have the approximate expressions given above [ffl]. Specifically, we have 
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Theorem 4.2. The Poincare map P has infinitely many fixed points labeled by 
I (I > lo): 



to — to,h x — &x 7 Q° — Qz\ z 2 — Z 2 J, 



where as I — > +oo, 



'(),/ 



-2.1 



i[/7r-^i] +o(l), 



} (i,0 



o(l), 



which x( 0,l \ Q^ ''' and z^ 1 ' 1 ^ are given in ^.6^,(^. r i),(4-i )■ 



5. Existence of Chaos 

One can construct Smale horseshoes in the neighborhoods of the fixed points 
of P. 

Definition 5.1. For sufficiently large natural number I, we define slab Si in 
S as follows: 



e s 



7 ? exp{-7i(i 0i2 (; + i) - ^r)} < 



z°(Q) < ?7exp{-7i(i , 2 ; - ^r)}, 
\x°(Q)\ < ?7exp{-ia t , 2 ;}, 
\4(Po{Q))\ <7?cxp{-ia<o,2;}, 
IIQ^oCQ))!! <*7exp{-ioi ,ai}}, 

where the notations x°(Q), 5|(Po (Q)), etc. denote the x° coordinate of the point 
Q, the z\ coordinate of the point Pq(Q), etc.. 



Si is defined so that it includes two fixed points pf and p^ of P (Theorem 4.2). 
Si, Pq(Si), and £P 1 (S';) are illustrated in Figure g, where £ is defined in (3.2). 
{ejo, e 5 o, e 5 o, e^o} denotes the unit vectors along (i°, zj 1 , z° , Q°)-directions in So, 
{e £ i , €yi , e 5 i , e^i } denotes the unit vectors along (x , y , z\, Q ^-directions in Si, 
and under the linear map £, {e^.i , egi , eji , e^i } are mapped into {fji , Syi , £ 2 i , £q! }. 
Let eg be the unit angular vector at Po(pf) of the polar coordinate frame on the 
(x 1 , y 1 )-plane. Let Eg = Leg, and we assume that 



(A3). Span<j e%o,eQ ,Eg,£ s i } = S 



Let Si l0 - — u o Si where the symmetry o~ is defined in (3.1). We need to define a 
larger slab Si such that Si U Si t „ C Si- 
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Figure 2. An illustration of Si, P<$(Si), and CP^(Si). 



Figure 3. (a) shows one of the homoclinic orbits, and (b) shows 
the blow-up of the neighborhood of the saddle Q e . 



Definition 5.2. The larger slab Si is defined as 



Si = Qes 



??exp{-7i(t 0! 2(( + i) - 2^)} < 



z°(Q) < »7exp{-7i(to,2; - ^)}, 
\x°(Q) - x°\ < rjexp{-^a t ,2i}, 
\4(Po(Q))\ < 14*1 +VeM~\a h, 2l }, 

IIQ^oW))!! <v<xp{-\ato, 2l }\, 
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where z\ „ is the z^-coordinate of Q\. 

Definition 5.3. In the coordinate system {x°, z®, z®, Q }, the stable boundary 
of Si, denoted by d s Si, is defined to be the boundary of Si along (x°, Q°)-directions, 
and the unstable boundary of Si, denoted by d u Si, is defined to be the boundary of 
Si along (z®, ^-directions. A stable slice V in Si is a subset of Si, defined as the 
region swept out through homeomorphically moving and deforming d s Si in such a 
way that the part 

d s Si n dJi 

of d s Si only moves and deforms inside d u Si. The new boundary obtained through 
such moving and deforming of d s Si is called the stable boundary of V, which is 
denoted by d s V . The rest of the boundary of V is called its unstable boundary, 
which is denoted by d u V. An unstable slice of Si, denoted by H, is defined similarly. 

As shown in [[jj], under the assumption (A3), when I is sufficiently large, P(«Sj) 
and P(Sz ](T ) intersect Si into four disjoint stable slices {Vi,^} and {V-i,V-2\ in 
Si. Vjf's (j = 1, 2, —1, —2) do not intersect d a Sf, moreover, 

(5.1) d s V t C P{d s Si), (t = 1, 2); 8 S V C P{d s S ha ), (i = -1, -2). 
Let 

(5.2) H j = p-\V j ), (i = 1,2,-1,-2), 

where and for the rest of this article, P" 1 denotes preimage of P. Then Hj 
(j = 1,2,-1,-2) are unstable slices. More importantly, the Conley-Moser con- 
ditions are satisfied as shown in pj]. Specifically, Conley-Moser conditions are: 



Conley-Moser condition (i) 



Vj = P(H 3 ), 

d s Vj = P[d B Hj), (j = 1,2, -1,-2) 
d u Vj =P(d u Hj). 



Conley-Moser condition (ii): There exists a constant < v < 1, such that for 
any stable slice V C Vj (j = 1, 2, —1, —2), the diameter decay relation 

d(V) < ud{V) 

holds, where d(-) denotes the diameter [§, and V = P(Vniffc), (fc = 1, 2, -1, -2); 
for any unstable slice H C Hj (j = 1, 2, —1, —2), the diameter decay relation 

d{H) < vd{H) 

holds, where H = P'^HnVk), (k = 1,2,-1,-2). 

The Conley-Moser conditions are sufficient conditions for establishing the topo- 
logical conjugacy between the Poincare map P restricted to a Cantor set in So, and 
the shift automorphism on symbols. 

Let W be a set which consists of elements of the doubly infinite sequence form: 

a = (• ■ -a_2a-iao, ai«2 ■ ••)) 

where S {1,2,-1,-2}; k e Z. We introduce a topology in W by taking as 
neighborhood basis of 

a = (• ■ •a_ 2 a_ 1 a ,a 1 a 2 • ••), 
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the set 



a k = a* k (\k\ < j) 



for j = 1, 2, • • •. This makes W a topological space. The shift automorphism \ is 
defined on W by 

X : W^W, 

Va S W, y(a) = where 6 fe = a fe+ i. 

The shift automorphism x exhibits sensitive dependence on initial conditions, which 
is a hallmark of chaos. 
Let 

a = (• • -0-20-100, ai«2 ■ ••)> 
be any element of W. Define inductively for k > 2 the stable slices 

K o_i — P(H a _ 1 ) r\ H ao , 

Va a-i...a-k = P\Ya-i...a-k) H -ff ao ■ 

By Conley-Moser condition (ii), 

d(V r O0O _ 1 ...a_J < ^l<i(V r „oa- 1 ...a_ (t _ 1) ) < ••• < d{V aoa ^) . 

Then, 

OO 

fc=i 

defines a 2 dimensional continuous surface in So; moreover, 

(5.3) 0V(a) C d u Si. 
Similarly, define inductively for k > 1 the unstable slices 

H aoai = P 1 {Ha 1 n T4 ), 

Ha a 1 ...a k = -P (Ha 1 ...a k H V^ ). 

By Conley-Moser condition (ii), 

d(H aoai ... ak ) < v 2 d(H aoai ... ak _ 1 ) < ... < v%d{H ao ). 

Then, 

OO 

#( a ) = P) H aoai ... ak 

k=0 

defines a codimension 2 continuous surface in Eo; moreover, 

(5.4) flff(o) C a s 5 ; . 
By (5.3|;5.4) and dimension count, 

V(a) n H{a) ^ 

consists of points. Let 

p e V{a) n H(a) 
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be any point in the intersection set. Now we define the mapping 

: W^Si, 



By the above construction, 



That is, 



(j)(a) = p. 
P{P) = Hx(a))- 

P o (f) — <p o Y. 



Let 



A = cf>(W), 

then A is a compact Cantor subset of Si , and invariant under the Poincare map P. 
Moreover, with the topology inherited from Si for A, cf> is a homeomorphism from 
W to A. Thus we have the theorem. 

Theorem 5.4 (Horseshoe Theorem). Under the generic assumptions (Al)-(A3) 
for the perturbed nonlinear Schrddinger system (1.1), there exists a compact Cantor 
subset A of Si, A consists of points, and is invariant under P. P restricted to A, 
is topologically conjugate to the shift automorphism \ on four symbols 1,2, —1, —2. 
That is, there exists a homeomorphism 

6 : W 



such that the following diagram commutes: 

W — i 



A. 



A 



A 



6. Numerical Evidence for the Generic Assumptions 

6.1. Generic Assumption (Al). Figure || shows a numerical result of Mark 
Winograd. It indicates that the homoclinic orbits are indeed tangent to the (x, y)- 
plane as t — > +oo. Specifically, Figure || (a) shows one of the homoclinic orbits, 
and Figure |H (b) shows the blow-up of the neighborhood of the saddle Q e . 

6.2. Generic Assumptions (A2) and (A3). Numerical simulation of the 
generic assumptions (A2) and (A3) is planned for a future work. 
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